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Abstract 

In this paper we take an idea presented in recent paper by Carlen, Carvalho, 
Le Roux, Loss, and Villani ([5]) and push it one step forward to find an exact 
estimation on the entropy production. The new estimation essentially proves that 
Villani's conjecture is correct, or more precisely that a much worse bound to the 
entropy production is impossible in the general case. 
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1 Introduction 



In his 1956 paper on the Foundations of Kinetic Theory ([5]), Mark Kac proposed 
a probabilistic model describing a system of N one dimensional, randomly colliding 
particles. The description is given by Kac's Master Equation 

dib 

-£{ Vl ,...,v N ,t) = -N(I-Q)i/>(y u ...,v N ,t) (1.1) 



at 



where 



(v 1 ,...,v N ) = V / 4>(R i)j (-&)(vi,...,v N ))d& 



with 

R i,0) ( V l> ■ ■ ■ > V n) = (Vi, . . . Vi{"9), . . . . . . ,v N ) 

Vi("&) = Vi cos $ + Vj sin Vj ($) = —Vi sin ■& + cos •d . 

The function ip(vi, . . . , vn, t) is a probability distribution on the energy sphere and it is 
formally given by 

^(. )t)=e -iV(/-Q) Vo 

for some initial condition tpQ. In the same paper, Kac introduced the notion of chaotic 
sequences (although he did not call it that way) and showed that this notion is preserved 
under the time evolution. This property is now called Propagation of Chaos. Kac went 
further and showed in fact that single particle marginal of the evolved density is a 
solution of the model Boltzmann equation 

df If f 27r 

— (v,t) = — dio dd (f (v cos # + uj sin??, t) f (— v sin$ + to cos i9, t) — f(v, t)f(u>, t)) 
dt 2tt J r J 

and thus giving a cogent derivation of the spatially homogeneous Boltzmann equation. 
For a detailed review the reader may consult [3]. 

The equation (II. lh . or rather the operator, is a bounded self-adjoint operator in 
the space I? (§ N ~ 1 (\ /r N), da N ^j where do N is the normalized uniform measure on the 

sphere. It is fairly easy to see that the time evolution defined by (jl.ip is ergodic, i.e., 
the solution will approach the function ip = 1 as t — > oo. By the spectral theorem, the 
rate of approach to the constant function in the sense of L 2 distance is governed by the 
gap 

A N =mf{(v,N(I-Q) V ) : {<p, 1) = 0, (<p,<p) = 1} 



where the infimum is taken over all <p e L 2 [ S^-^ViV), dcr ■ Kac conjectured that 



lim inf Ajv > . 



2 



The conjecture was proved to be true by Janvresse in ([3]) and the exact value of A^v 
was computed by Carlen, Carvalho, and Loss in (|2j). 

The L 2 distance is rather unsatisfactory. For any reasonable density ip, in particular 
a chaotic one, it is easy to see that 

U(vi,...,v N ,0)\\ L ^ sN _ 1{VN)>daN > j > C N 

where C > 1 and hence it would take a time of order N to see a substantial decay of the 
I? . Clearly, this is not what one considers "approach to equilibrium". A more natural 
quantity to use is the entropy 

The crucial difference between the I? distance and the entropy lies in the extensivity 
of the entropy, namely that if V'w (v i, . . . , vjy,t) satisfies ip^ (v±, . . . , v n, t) ~ n^ 1 /(^j, t) 
in a weak sense, i.e., chaotic (referred by Kac as 'The Boltzmann Property') then 

H N {il; N ) « nJ f(v tt ) log (^^) dv = NH(f(v,t)\ 7 (v)) 

where j(v) is the normalized Gaussian. 

Differentiating the entropy of a solution to the Kac Model gives the time evolution 
equation: 

dH N (il) N ) 

— = (log ip N ,N(I - Q)tp N ) 

This, along with a known inequality by Csiszar, Kullback, Leibler and Pinsker and 
the enxtensivity property allows us to conclude that 

. . .,v N ,t)da N - da N \\ 2 Tot&l Variation < 2N e - r " t H(f(v, 0)| 7 («)) 

for 

. f (log(j; N ),N(I-Q)^ N ) 
1 at = mt TTTT~\ 

where the infimum is taken over all probability densities ipjy on E> N ~ 1 (\fN) which are 
symmetric in all their components. Tat is called the entropy production. 

The hope that there exists C > such that r^r > C was refuted in 2010 in an paper 
by Carlen, Carvalho, Le Roux, Loss, and Villani ([3]) where the authors managed to 
find a sequence of probability densities {0Ar}^r gN with 

r (log(0 j y),iV(/-Q)0 j v) 

hmsup = (1.2) 
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While this means that the time of convergence to equilibrium is not of logarithm type, 
an exact estimation on the entropy production might still give a better convergence rate 
than that of the original Kac model. 

The first step towards this goal was done in 2003 by Villain in (|6j) who proved that 



N 



> 



N - 1 

Villani conjectured that 

r " = ° ik 

which wouldn't bode well for the approach to equilibrium in the ergodic sense, but poses 
an interesting mathematical problem. 

The main result of this paper is to show that Villani's conjecture is essentially true. 
More precisely, we will show that 

Theorem. For any < f3 < i there exists a constant Cp depending only on (5 such that 
(See Theorem [T7] in Section @|. 

Both (ll.2p and (I1.3P are proved with the same idea: creating an N particle symmetric 
function F^ from a one particle function / 

*N (Vl, ■ ■ ■ ,V N ) 



ZnU, VN) 

where 



Z N (f,r)= [ Uf =1 f( Vi )dai 



i(r) 

and da^f is the uniform probability measure on S Ar_1 (r). The main difference between 
the two proofs lies in the fact that while in ([3]) / remains fixed, in our paper / changes 
with N via a parameter 5 = 5m- 

The paper is structured as follows: Section [2] reviews known results about the nor- 
malization function Z^(f, t). Section [3] is our main theoretical part of the paper, dealing 
with general properties that will allow us to give an asymptotic expression to the nor- 
malization function. Section [3] is where we prove our main result. Picking a function 
which is natural to the problem at hand and using the result of the previous sections 
along with some involved computation. Section [5] contains a few last remarks and the 
Appendix has some simple but very useful computation that we use throughout the 
entire paper. 

We'd like to conclude the introduction by thanking Michael Loss for his helpful 
remarks and discussions, making this paper possible. 
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2 The Function Z N (f,r) 

The key to the computation of the entropy production lies with the normalization func- 
tion Zpj(f, r). In this short section we'll find a simple probabilistic interpretation to it, 
along with a formula that will serve us in the following sections and the final computa- 
tion. This section is a short review of known results from ([3]). 

Lemma 1. Let f be a density function for the real valued random variable V . Then the 
density function of the random variable V 2 is given by 

Proof. For any function ip = (f(\x\) = ip{r) we find that 



Eip= / y(r) •(/(?■) + /(-r))dr 
Jo 

on the other hand 

Eip = J ip (\/t) h(t)dt = J cp{r) -2r-h (r 2 ) dr 

Since (p was arbitrary we find that 

2r • h (r 2 ) = f{r) + f(-r) 

and the result follows. □ 

Lemma 2. Let V\, . . . ,Vn be independent real valued random variables with identical 
density function f(v). Then the density function for Sn = Ei=i ^ s <?2i>en by sn{u) = 

^uf-lZ N (f,^). 

Proof. Similar to Lemma Q] for any (p = <p(r) we find that 

\ roo 
N I j_ / _/„\ioiV-l|.JV-l, 



%>=/ <p(r)\ f(vi) ■ ■ ■ f(v N )ds J ; )dr= <p(r)\S"- L \r"- L Z N (J,r)dr 

J0 yjS^-^r) J Jo 

on the other hand 

Eip = / f(\/x)sN(x)dx= / <p(r) ■ 2r ■ sjy (r 2 ) dr 
Jo Jo 

Since <p is arbitrary 

2r. SN (r 2 ) = \S N - 1 \r N - 1 Z N (f,r) 
which implies the result. □ 
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Corollary 3. (Expression for Z]\f(f,r)) Under the conditions of Lemma\^ 

2h* N {r) 



Z N (f, VP) 



|AT-l| r f-l 



where h N is the N-fold convolution of h, defined in Lemma[l\ 

Proof. This follows immediately from Lemma [21 Lemma [T] and a known probability 
fact. □ 

3 Central Limit Theorem 

In order for us to be able to compute the entropy production an asymptotic behavior 
for Zw(f,r) is needed. As seen in Section [2] the function Z/\r(/, r) is closely related to 
the A-fold convolution of the density function h{u) and as such we'll employ standard 
techniques to estimate it. The specific function we'll construct as a test function for 
the entropy production has the property that the Fourier transform of its one particle 
function splits the line into two natural domains: One where we can use analytic expan- 
sion, and one where the decay is dominated by exponential functions. The radius of the 
separating circle would depend on a parameter 5 =<5jv that we'll exploit later on to get 
the final conclusion. 

While this is the case arising in our specific construction, we believe that it's a 
natural way to view the problem. Even though we have yet to attempt any different 
test functions we think that similar situation would happen in a larger class of functions 
created from one particle function. As such, a generalization of our computation was 
made and is presented in this section. 

The reader should keep in mind the following intuition while reading this section: 
<?(£) represents the Fourier transform of the function h(u), connected to the one particle 
function via Lemma [TJ The first lemma of the section explores the domain outside the 
radius of analiticity while the second explores the domain where analytic expansion is 
possible. Lastly, the parameter^ is a function of N, going to zero as N goes to infinity. 

Lemma 4. Let gs(£) = gs N (£) be such that 

00 f or l£l > C( 5 \ds{0\ < 1 — ot(6), where a(5) > 0. 

(*») \9S(0\ < 1 for allt 

Then 

[ |sf(e)-7f(0|de 

J\e\x* 

— — i 

J\Z\>c6 7rc<5£| 7rcdS| 

where 7l (f) = e" 2 ^ • e~ 2 ^^l . 
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Proof. We have that 



sf(0—^(0|de 



\Z\>cS 
N-l 



\gs(0-Ji(0\ 



N-l 



k=0 



l«l> rf fc=0 

2V-1 



< 2 [ IgsiOf- 1 dt + 2 £ (1 - atf))"-*" 1 / e" 2 ^^ 



'\t\><* k=1 
Using Lemma [18] and [19] in the Appendix we find that 



N-l 

E 

k=kg 



151 > c5 k=k 



N-l /=- 

Z7T • e 2 



4/br 2 £2 



< 



2tTC<5£ 2 ; 



-2fc 7r 2 c 2 5 2 E 2 



Hence 



/ K(e)-7f(0|^ 

7|^|> C 5 



< 2 



ICI>c<5 



l^(Or _1 ^ + 2(l-a(5)) 



N-l 



+ 2 E 

fc=[f]+1 ^l«l>c<5 



N-[f)-lJ2 

=1 



fc=i 

2fc7T 2 £ 2 E 2 



-2fe7T 2 ^ 2 S 2 



ICI>c<5 



< 2 



m>cs 



i gi «)r-'< + (1 ' a( gr' + ' 



7TC(5S| 



-(l+JV)7r 2 c 2 5 2 S 2 



Lemma 5. Let g${C) = gs N {Cj he such that 



□ 



(t) there exist M , Mi, M 2 > such that sup| ? | <c5 \g s (£) - 71(f) | < (fR + + M 2 ) |£| 3 . 
(ii) for c5 1+l3 < |£| < c<5 < 1 — w/iere 0^(5) > 0. 

(m) 1^(01 < 1 for all 
Then, 

I \ g » {i) -,nm<- cH2(Ma+ T +M ^ 

J\t;\<c5 2 
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cHVN (Mo + MiS + M 2 S 2 ) (1 - a^d))^ 1 cH 1 ^ (M + Mrf + M 2 £ 2 ) e -* 2 (N-i)<?S*+W s 

+ 2c 3 (M + Mi 5 + M 2 S 2 ) VN5 1+3/3 
where 71(f) = e" 2 ^ • e ~ 27r2 « 2s i. 

Remark 6. The coefficients Mo, Mi and M2 play a major role in the estimation. Notice 
that we can get a better result if have that Mo = and an even better result if both Mo 
and Mi are zero. 

Proof. Similar to Lemma U] we find that 

r N ^ r 

/ bf(o-7f(6K<E / ^s(o - 7i(oi i^or - *- 1 iTiCOi** 

J\t\<cS k=Q J\(\<c6 

2 y c «5i+'5<i5i< C 5 v d d / 

^ yi€i<c5i+/9 v ^ 2 <* / 

We have that 

E/ (§ + ir + M2 ) i^i 3 ^(or*- 1 bi(e)i* 
y c <5i+/3<i5i< C 5 v ^ / 

iV-l „ 

< c 3 5 (M + Afi* + M 2 5 2 ) Y (1 - a^O*))*"*" 1 / e" 2 ^ 2 ^ 

[f] . 

< c 3 5 (M + M X S + M 2 5 2 ) (1 - a^))^ -1 V / e" 2 ^ 2 ^^ 

N-l r 

+c i 5(M + M 1 5 + M 2 5 2 ) Y / e - 2fc ^ 2s2 ^ 
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< c 3 S (M + Mi 5 + M 2 8 2 ) (1 - a^S))^" 1 ^ 



k=l 

N-l 



+c 3 5{M + M 1 5 + M 2 5 2 ) y If e ~ 2k ^ 2 sdZ- [ 

, 7Z7\ . \J\£\<c6 J\t 



m<cS J\(\<oS^+P 



fc=lf 1+1 



L Z J -I 

< c 3 <5 (M + Mi 5 + M 2 5 2 ) (1 - a^))^" 1 —== 

fc=i J2nY? s k 

+c 3 5(M + M 1 5 + M 2 5 2 ) ^ 



< 



2tt£ 2 



1 




"TV" 












~2 



| c 3 <5(Mq + MiJ + M 2 J 2 ) g 1 e -2.W+^l_ e -i.Wi;l 



„ c 3 5^N (M + Mi 5 + M 2 <5 2 ) (1 - a p {5))^~ 1 



c 3 <5(M + Mi5 + M 2 <5 2 ) ^ e -^fcc»^E» 



2.S 2 Vl-e-^ 2 ^ 



^ c 3 5y/N (M + Mi (5 + M 2 5 2 ) (1 - (<$))"* 



-i 



^(Mq + M^ + M^ 2 ) ^ e -2^fcc^+^sg 



l+i 



< c 3 8VN (Mq + Mi<5 + M 2 5 2 ) (1 - c^))^' 1 | c 3 5 1 -? (M + Mi*5 + M 2 5 2 ) e -^(N-i)^ 2 + 2 ^ 2 



vrS 2 2vr C( 5S 2 • Vl - e' 2 - 2Nc2s2 ^ 
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Next we find that 
N-l 



i j\z\<c&^+p 



^ + f + M 2 )\e\9s(0\ N - l '- l bm\"^ 



N-l 



< c 3 (M + M x 5 + M 2 5 2 ) 5 1+3f3 • V / e^ 2 ^ 2 ^ 

i|5|<c5i+/3 

AT-1 

< c 3 (M + MiJ + M 2 8 2 ) 8 l+w ■ 



1-e 



J -4fc7r 2 c 2 <5 2 + 2 / 3 S2 



k=l 



< 



c 3 (M + Mi (5 + M 2 5 2 ) 5 1+3 P 



N-l 



< 



■y- 

2c 3 (M + M 1( 5 + M 2 5 2 ) ^/N5 l+ ^ 



Which completes the proof. 



□ 



Theorem 7. Let h§(x) = hs N (x) be a function such that gs(£.) = hs(£,) satisfies 

(i) for |£| > c5 N \gs N {Q\ < 1 - a(6jv), w/iere a((5jv) > 

(ii) there exist M , M u M 2 >0 such that su P | 5 | <C(5iV \gs N (£) - 71 (01 < (if + fr + M s) 
(m) /or co"^ < |£| < c5at 1^(01 < 1 - <*p(8 N ) where ap(5 N ) > 

( vi ) \98 N {Q\ < 1 M all £ 
and if 

5n, cv(5n) andap{5j<j) are domianted by powers of N 



a(5]\r)N — > 00 
(xb(Sn)N — > 00 

' N->oo 

WN-> 00 

JV N— >oo 

S ]+^ N 

JV JV->00 

o"^ 1 ^ T<s N is bounded 



(3.1) 



then 



sup 



N 



'NE Sn V2vr 
where h* s N {x) is the N-fold convolution and e(N) 



< 



e(N) 



NZ S , 



JV-s>oo 
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(x-N) 2 
2NS'4 

1 e o 



Proof. It is easy to check that v /-^ s • e ^ 2 - (£) = 7 1 v (£) 
Using Lemma 0] and [5] we find that 



sup 

X 



hf(x) 



0-AQ z 

e 2JVE * 



iVS 5 V2vr 



< 



jv 



(0-7f (0|« 



^|e|<c<5 ->|£|>c5 

1 / c 4 VN5 1 +W6% ( - 1 -ft'E5 (M Q + M x 8 + M 2 8 2 ) 



NZ 5 \ 2 

c 3 5iV (M + Mi5 + M 2 S 2 ) (1 - (PVNS 1 -? (Mq + + M 2 <5 2 ) e -7r 2 (iv-i)c 2 .5 2 + 2 ^g 



7T 



+ ^ ' ^ + 2^2,5 / I^COI^- 1 

V27T J|?|>c<5 



_ (1+Ar)7r 2 c2(5 2 E 2 



A/ 



+2(1 -a{5))^- 1 ~= + 



27TC<5£,5 7rC<5£,5 

Conditions (13.11) insure the desired conclusion. 



□ 



Remark 8. A careful look at the proof of TheoremOshows that for a fixed j if lini7v->oo V~N — j^s N f\£\> c 8 N \d$N 
and conditions (|3.ip are satisfied (with the obvious change) then 



sup 

X 



CM 



(x-N+j) 2 
2(N-j)S 2 



< 



6,(JV) 



where eJN) — > 0. 

J iV-^oo 



4 Entropy Production and Villani's Conjecture 

In this section we'll find an exact estimation for the entropy production. The idea 
behind this estimation is to use superposition of stationary solutions for the Boltzmann 

equation: the Maxwellian densities M a (v) = ^ 2 == . This idea was exploited by Carlen, 
Carvalho, Le Roux, Loss, and Villani (|3j) and Bobylev and Cercignani ([!]) before them. 
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The basic one particle function would be 

fsM = fs(v) = 5Mi(v) + (1 - S)Mi(v) 

Id 2(1 — o) 

This function has the property that both its parts have the same energy 
[ 8MUv)dv= [ (\-5)M, (v)dv = \ 

JR 25 JR Al-<5) I 

while as 5 gets smaller the number of particles represented by SMj_(v) is far smaller 
than those represented by (1 — S)M i (v). The fact that we have a small number of 

very energetic particles and a large number of very stable particles trying to equilibrate 
will cause slow decay into equilibrium. That physical intuition is indeed true as would 
be seen shortly. 

Lemma 9. Let h s (u) = Mv^l±MzVg = My^) then 

(*) f °° h s (u)du = 1 
(ii) uhs(u)du = 1 

(Hi) £ 2 = J °° u 2 h s (u)du - (/ °° uh s (u)du) 2 = iS{ ^_ s) - 1 



l-<5 



Proof, (i) — (Hi) follow immediately from the fact that J °° u m hs(u)du = J m x 2m f$(x)dx 
and the fact that 

/ M a (u)du = 1, / u 2 M a (u)du = a, u 4 M a (u)<iu = 3a 2 
Jr Jr Jr 

We're only left with proving (iv). 
It is easy to check that 

± f M a (u) • e- 2 ^ 2 du = - ~ 2 7 m - / M a (n)-e- 2 ^ 2 ^ 
The initial value problem JfV^O = i+^^ fiO^ £ £ ^> ^(0) = 1 nas the unique solution 

m 1 



a/1 + 47rla| 

Thus, the result follows from the definition of and the fact that 



M0= r h & (u)e- 2 ^ u du= [ f s (u)e- 2 ^ u2 du 

JO JM 



□ 
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Lemma 10. Let gs(£) = hs(£) where 5 < \ then 



2 

i/4 
5 



p\ (S) where pl {^ 



(i) for\Z\>± \ gs ^)\<l-s{l- 

(ii) there exist M , M 1: M 2 > such that sup m<A _ \g s (£) - 71 (£)| < + ^ + M 2 ) |£| 3 . 



5-s>0 



(Hi) for < \£\ < £- \g s (0\ < 1 



An I s I 47T 
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+ />2(5) tf/iere # 



(«) /or a fixed j Jj^jL 1^(01 5 d £ < 



iV-j-l 



<$->-o 

N-j-l 



+ 



n(N-j) 



Proof, (i) For |£| > £ 
l<w(0l < - 



1-5 



4 /l + 4^i 



< 



1 + 



4?r 2 g 2 



+ 



1-5 



4/5 4 /1 + 



= { -s + (i-s) 1- 



5 2 



16(1 - 5) 2 



+ . . . 



1-5 1 



£ 2 

4(l-5) 2 



5 +«M 



where ^ — > 0. 



(m) Using the expansions for and e x we find that for |£| < ^ 



IM0-7i(0l<KI 

3vT 3 



/Stt 3 
I 5 2 



(¥) 



+ 



8vr 3 



5(1 - 5) 
+4vr 5 



4^ 3 + 2vr 4 
3 



45(1 - 5) 

2 

2 



(1 - 5) 2 

37T 4 



5(1 - 5) 



1^1-4^1^1 



^45(1-5) 
+87r 3 |^(-27rif)| + 8vr 6 



1) 1^ + 4^ 
3 

45(1 - 5) ~ 



45(1 - 5) 



1 itf 



1 |fl s |^(-2**Efc 2 )| 



where </>(a;) is analytic in \x\ < \ an( i VK X ) is an entire function. Denoting = 
\<t>{x)\ and = sup| x |<i we find that 



|^(0-7i(0l < 
(in) For |d > ^ 

\98(0\ < 



/8vr 3 57tt 3 o/ 



1-5 



1 + ^ 



1 + 



47T 2 g 2 



< 



+ 



1-5 



52+2,3 
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where j-t 0. 

(iv) This is a general property of the Fourier transform of a density function. 

(v) 

N-j-l 

JV-i-l 

H , dx 



l— \ \ N-j-l 



\ 



7T 7r(iV - j - 3) 



Remark 11. Note that in our case 



•/|£|>c<% 



N-j-l 



< ' ; + 



□ 



2n6 VN6-yJl-i# 

so as long as the conditions in (|3,ip are satisfied we have that e,-(iV) defined in Remark 
E] would satisfy e,-(JV) — ► 0. 

Theorem 12. Lei /^(u) = = 5Mj_(v) + (1 - <5)M_i (u) suc/i i/W 

25 2(1 — 5) 

5 at is domianted by powers of N 

6 N • Ar ^ 00 00 (4.1) 
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then for a fixed j 



z N-j {fs N ,Vu} 



/ (n-iV+j) 2 

e d N 



N—i i 



\ 



2tt 



+ A j (iV-j» 



V 



/ 



where sup ugK |Aj(iV — j, u)\ < €j(N) and limjv-»oo ej(N) = 0. 

Proof. This is immediate from Lemma [31 El EH Theorem [7] and Remark [TT] 



□ 



We're now ready to compute the entropy production. We'll start by estimating its 
denominator and numerator. 



Lemma 13. Let Fn (v±, . . . ,vjf) = where 5^ satisfies conditions l\4-l\ )- Then 



lim 

7V-s>oo 



h"-HVN) f n log F N da N i g 2 



N 



Proof. Using the symmetry of the problem, Lemma [22] from the Appendix, Theorem [T2l 
and Stirling's formula we find that 



/ F N logF N da N = 1 -fl f fiiLifsivi)) log f 5 (v k )da N -log Z N (f s , 

Jsm-i(Vn) Z N (f s ,VN) ^Jsn-hVn) 



X\S N - 2 \ , -2=2 Z N-1 (f6,V N - v i) r- 

AWlog/iW iV-^ 2 ,V J d Vl -logZ N (f s ,y/N) 



N 



/i(«i)log/i(t;i) -X[-y/N^(vi) 



■ e {N ~ 1)s 'i + V2^Xi (N — 1, N — v\) dvi 



log (V2 + O f-^Lj ) (l + ^A (iV, JV) 



iV 1 
y (log 2^ + 1) ---log 



Since < fg < 1 we have that 



fsivt) log fg^-x^^vx)- e (Ar - 1)E i + v / 2^A 1 (iV-l,iV-?; 1 2 ) 
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< (l + V2^ei(Nfj (-fsM log f s (vi)) 
<(l + ^ei(If))(-8M±(v 1 )]ag(6Mj L (vx))-Q.-8)M i («i) log ((1 - S)M i 

\ / \ 25 \ 25 / 2(1 — 5) \ 2(1- 

= 35(«i) 

It is easy to check that 



9s N (v) — > -Mi (u) log Mi (w) 

AT— S>0 2 2 



and 

log 7T 1 



/ gs K (v)dv — > - Mi(v)logMi(v)dv 



2 + 2 



Since 



( i(l-vl) 2 S N (l-S N ) 
e (iv-D(3-^(i-^)) + V2^Ai (N-1,N- v\) 

— > Mi (vMogMiivi) 

N-^OO 2 2 

we conclude that 

h*-i(VN) Fn 1o S F Nda N . , , 1 log2vr log 2 



/ Mi(iri)logAfi(«i)dt;i + - + 

JR 2 2 ^ 



iV at-xx, y R 5 V x/ ° 5 V ±y 2 2 2 

due to the generalized dominated convergence theorem. □ 

Lemma 14. Let Fn (v i, . . . , vn) = ^~\f ^/W) w ^ ere ^ N sa ^ s fi es conditions f//. Then 
there exists a constant Ct yve -8 depending only on the behavior of 5jy such that 

(log F N ,N (I -Q)F N ) 



N 



< C type _s (S N log 5 N ) 



Proof. Similar to Lemma [T3l by using the symmetry of the problem, Lemma [221 from the 
Appendix, Theorem [12] and Stirling's formula we find that 

(log F N ,N (I -Q)F N ) 

1 N f 

Z N (f s ,VN)(N -1)tt iHJ&x-i(VN) 



\i<j j ° 



vn) - f® N (RiM K ■ ■ ■ , vn))) d-d da N 



16 



if i and j are different than k the integral is zero and so 

N 



(log F N ,N (I -Q)F N ) = — * / tog 

^(/5,VA^)(7V-l)7r^^7s JV - 1 (ViV) 

^ (/ 0JV K • • • , VN) - /®" (flfcjW K • • • , 

AT ^ 







^ / (-log (/«M<?))/«M0)) - fs(vi)fs(v 2 )) {nlsfsivi)) da N 

f'ZTT r 

IT JO Jv\+vl<N 



ZnUs, Vn)tt Jo Js n - 1 (Vn) 



2 2 ,*=±Zn-2 (fs.y/N-v'i-vi) 



2= r d{> I (-log (h(vimh{v 2 m - fs{vi)fs(v 2 )) 

l_lio Jvl+vl<N 



e (JV - 2)E ^ + ^A 2 (iV - 2.N -v\- vl) 
== dv\dv 2 

Using rotational symmetry and symmetry in the variables we find that 

(log F N ,N (I -Q)F N ) 

= ^= C 1 * dti [ (logf S (vimfs(v2m - logf S (vi)fs(v2)) 

4vr Jl - jj J o Jv\+vl<N 

_(2-v{-vl) 

(/*M*))/*(i*W) - /,(i*)/,(t*)) 7=^— * — 2j -d Vl dv 2 

1 + V^7rAo(iV, iV ) 



< ^= r * d$ [ (log fsiv^fsiviW-lozhin .)f 6 (n)) 

4vr^l -%Jo J* 2 



1 + v^7i"Ao(iv, iv J 



AT (1 + ^e 2 (iV)) 



irJl-%(l + y/to\o(N,N)) 



f 2W d$ [ (-log f s ( Vl )) (friv^fs^)) - fs( Vl )f S (v 2 )) d Vl dv 2 

Jl 2 
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Since M (ui(i?))Af (u 2 (i?)) = M a (vi)M a (v 2 ) we see that 



fstyiWfsfaWyfsMfsfa) = 6(1-5) Mi (di(i?))M i (^(,?)) - Mjl(i*)M _i fo) 

\ 2(5 2(1 — (5) 26 2(1 — (5) 

+5(1-5) (Mi (v 2 (-d))M i -Mi (v 2 )M i 

V 25 2(1-5) 24 2(1-4) / 



< 5(1-5) MlfafflW 1 fa(t?))+Mlfa(j?))M 1 

V 24 2(1-4) 24 2(1-4) 

and along with 

- log/*(ui) < - log (*Mj. (vi)) < + ^ + <5 + U2 2 (??) ) 

we conclude that 



(log F N ,N (I -Q)F, 



N 

TrT (1+VSAo(JVAr)) ; Ra v 2 2 ' 

< + f- - 1?£^ ! U ( -i,o gi ) 

'l _ 2. (1 + V2?A (JV, JV)) V2 21ogi 21ogi 21ogi; 

The result follows. □ 

Theorem 15. Let Fn (v\, . . . , vn) = ^{f ^/n) w ^ iere °~N satisfies conditions J^. 
Then there exists a constant Ct ype -5 and an integer N type _s depending only on the be- 
havior of 5n such that for every N > Nt ype -S 

(log F N ,N (I -Q)F N ) . i n _r \ 

< ^type-5 (-ON log 5n) 



Proof. This follows immediately from Lemma [131 and [141 □ 

Theorem 16. Let F N (v ly . . . ,v N ) = -^J^T where 5 N = jp^w and < P < h 
Then there exists a constant Cp and an integer Nr depending only on f3 such that for 
every N > Ng 

(log F N , N(I - Q)F N ) < C p log N 



I s n -i (VW) F N log F n da N ~ 



Proof. This follows immediately from Theorem 1151 and the fact that 5n = N i-2p satisfies 
conditions (|4.ip . 

From this we conclude our main result: □ 
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Theorem 17. For any < (3 < | there exists a constant Cp depending only on /3 such 
that 

CplogN 
Tn - N i-2P 

5 Final Remarks 

One question we might ask ourselves is: Can we modify the given proof to get the exact 
value in Villani's conjecture? Looking at the proof we notice that the result we obtained 
has very tight conditions in terms of j3. We needed 5^~ 2 ^ N to diverge to infinity and 
S^N to go to zero. This doesn't leave much room for variations. This leads us to 
believe that the family of functions constructed here would not be helpful to prove the 
exact version of Villani's conjecture. Something more clever must be done. 

Another question we don't know the answer to is the fourth moment question. Both 
in this paper and in (|4J) the family of functions constructed has an unbounded fourth 
moment. Would restricting the fourth moment lead to a lower bound on the entropy 
production? 

Lastly, can our computation be generalized to a more difficult interaction than Kac's 
model? Can we try and use the same idea in a different models of the Boltzmann 
equation? 

While we don't know the answers to the proposed questions we hope that this paper 
shed some light on the entropy production problem and that at least some of the above 
questions would seem more solvable after reading it. 



A Helpful Computations 

The appendix consists of Lemmas that are vital for the computations needed in our 
paper, and are used extensively in Sections [3] and 31 

Lemma 18. (Gaussian Integral Estimation) 



2tT I ay 2 f a 2 x 2 \>2lT I 5~ T 

1 - e~— < I e — ~dx < VI - er a r ? 

\x\<7] 



a 2 x 2 V27T ' e 2 

e 2 dx < 



\x\>7] 



Proof. We have 



f _a 2 jP L , f f a 2 (x 2 + y 2 ) / f f a 2 ( X 2 +y 2 ) 

/ e 2 dx = a e 2 dxdy < \ e 2 dxdy 

!\x\<r) V J J\x\,\y\<7] V J Jx 2 +y 2 <2ri 2 
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Lemma 19. (Special Sums Evaluation) 



g 2 fc , _ a fe() 

e " 2 V2vr • e 2 



< 



K = K0 + 1 



Proof. We have that 



> — 7=~~ - / — ^cfe = - / _ e 2 dy < - 

k=ko+l Jko V x y=a\fx a J aVh) dJaVh) 



1 — a K o 

If _i£ , v^-e — 2- 

= — / e 2 ay < 

a J\v\ 



'\y\>av% a 

Similarly 



E 4< ^ = 2( v ^-v / ^)<2v / ^ 



The next set of Lemmas refer to integration over the sphere § 1 (r 
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Lemma 20. (Integration on the Sphere I) Let f (v\, . . . , fjv) be a continuous function 
on M. N then 



/ ffo? = V / ^ . J -dv x . . . dv N . 



Proof. Standard in any Differential Geometry course. □ 
Corollary 21. (Integration on the Sphere with the Uniform Probability Measure) 



dv-i 



Lemma 22. (Integration on the Sphere II) Let f (y±, . . . ,Vj) and g . . . , Ujv) be 

continuous functions on W and W N ~ J respectfully. Then 



f (vi, ...,Vj)-g (v j+1 , ...,v N ) da? 



JV-j-2 

3 \ 2 

V? 

i=l 



/ , s gda N A- : I dv\ . . . dv 



Proof. Using Corollary 1211 we find that 

/ (vi, ...,Vj)-g (v j+1 , ...,vn) da? 



={+,-} 



f(vi,...,Vj)-g (v j+1 , v N -i,eyJr 2 -^Ji vf 



S N-l\ r N-2 J N _ x I „ 



(foi . . . dvN- 



|§JV-l| r Ar-2 / 



| S JV-i-l| 
| giV -l| r iV-2 



N-j-2 



3 \ 2 
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Is 



§N-j-l 



-?) v^ 7 




cfoi . . . (if j 



□ 
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